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In this paper we study the interaction of a cosmic string with a (non-rotating) black 
hole. This problem is interesting for many reasons. First of all, it is a rare example of 
gravitational interaction between two extended relativistic objects that allows complete 
analysis. When a cosmic string comes close to a black hole and is captured by it, one 
might expect generation of strong gravitational radiation, which makes the problem 
potentially interesting for astrophysical applications. Besides this, the study of the 
interaction of a cosmic string with a black hole has demonstrated some unexpected, 
physically interesting possibilities, such as the creation of 2D black holes [lj. 

For astrophysically interesting cases, the dimensionless parameter fi* = Gfi/c 2 , 
where fi is the string tension, is small. For example, for strings formed during a GUT 
phase transition, fi* pa 1CU 6 . For this reason, in the leading order one can neglect 
gravitational backreaction effects and consider the string as a test object propagating in 
a given gravitational background. The worldsheet of such a string is a minimal surface 
which gives an extremum to the Nambu-Goto action ||, |3|. 

Stationary configurations of a cosmic string in the field of a rotating charged black 
hole are described by worldsheets which are tangent to the Killing vector generating 
the time shift. Because of this symmetry, the problem is reduced to finding a geodesic 
in a three-dimensional space; the latter can be solved by separation of variables ||. In 
the absence of such a symmetry, non-linear partial differential equations describing the 
string motion in the gravitational field of the black hole do not allow exact solutions. 
Under these conditions one either needs approximation schemes or numerical solutions. 
In the present paper we assume that the string is initially far from the non-rotating 
black hole, straight, and moving with velocity v in the direction of the black hole. The 
motion of the string is affected by the gravitational field of the black hole. For a given 
velocity v, the outcome of the interaction depends on the impact parameter b. If the 
impact parameter is large enough, the string escapes the gravitational attraction of the 
black hole and is inelastically scattered. For smaller values of b the string is captured 
by the black hole. 
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The scattering of a cosmic string for very large impact parameters was studied in 
[7J. Such a string passes far from the black hole and its motion is modified only slightly. 
Under these conditions one can consider the gravitational interaction as a perturbation 
and treat it by expanding the solution in powers of the Newtonian gravitational potential 
(p = GM/(Rc 2 ). The numerical study of string scattering by a black hole was initiated 
by Moss and Lonsdale || . In Reference |§ we used numerical calculations to obtain the 
dependence of the critical capture impact parameter on velocity v, and demonstrated 
that the results of Moss and Lonsdale do not correctly reproduce the behaviour of the 
critical impact parameter at high velocities. 

In order to obtain more detailed information on the characteristics of scattering and 
capture of a cosmic string by a black hole we performed further numerical calculations. 
In regimes that allow analytical study, we compared numerical and analytical results. 
The agreement of these results was used as an additional test of our calculations. We 
also analysed the role of tension in the string by comparing the relativistic string results 
with those for a string without pressure (a so-called "dust model"). These results are 
discussed in the following sections. In Section 2 we formulate the scattering problem for 
the straight string. Section 3 outlines the numerical method that was used and discusses 
the scattering at shallow impact parameters. It also contains comparison of these results 
with the results obtained in the weak-field approximation. Section 4 compares numerical 
and approximate analytical results for a string moving with ultra-relativistic velocity. It 
also discusses coil formation for the scattering in this regime. (We make the distinction 
here between coils, which are loop-like configurations in a straight string, and the term 
loop, as found in the literature, which designates closed, bounded string configurations. 
We thank an anonymous referee for suggesting this clarification.) Section 5 discusses 
near-critical scattering, where the string worldsheet is partially wrapped around the 
black hole. Section 6 discusses the role played by the tension of the string and makes a 
comparison with the model of a dust string. Section 7 discusses aspects of the problem 
connected with string capture. 
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2. Scattering Problem for a Straight String 

2.1. Equations of motion 

The two-dimensional worldsheet representing a moving string is parametrized as X^(( A ); 
X^ (/i = 0, 1, 2, 3) are the spacetime coordinates, and ( A (A = 0, 3) are the worldsheet 
coordinates (° = t, ( 3 = a. The equations of motion of the string follow from the 
Nambu-Goto action [[Kj which we write in the Polyakov form [|J 

I[X,h] = -- [ dTdaV^hh AB G AB . (2.1) 
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We use units in which G = c = 1, and the sign conventions of fil l. In ( |2.1|) Hab is the 



internal metric with determinant h, and Gab is the induced metric on the worldsheet, 

Gab = 9^~qFa~q(b = 9^ x> X x ^b , ( 2 - 2 ) 
where g^ is the spacetime metric. 

The equations obtained by varying action (|2.1| ) with respect to X^ and Hab are of 
the form 

OX" + h AB T£ g X'XX f} B = (2.3) 



and 



GAB-\h A Bh CD GcD = 0, (2.4) 
where □ = (— h)~ 1 / 2 dA [(— h)~ l l 2 h AB dB\ ■ Equations (|2.3| ) are the dynamical equations 



for string motion, while equations (2.4) are constraints. 



2.2. Straight string motion in flat spacetime 

Consider first the motion of a straight cosmic string in the absence of an external 
gravitational field, that is when = rj^, where rj^ is the flat spacetime metric. In 
Cartesian coordinates (T, X, Y, Z), T]^ = diag(— 1, 1, 1, 1) and T^g = 0, and it is easy to 
verify that, taking h A B = Vab = diag(— 1, 1), 

X" = X"(r, a) = (cosh(/3) r, sinh(/3) r + X , Y , a) (2.5) 

is a solution of equations fl2.3p and (|2.4|). This solution describes a straight string 



oriented along the Z-axis which moves in the X-direction with constant velocity 
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v = tanh/3. Initially, at r = 0, the string is found at X^(0,a) = (0, X ,Y , a), where 
Y = b is the impact parameter. For definiteness we choose Y > and X < 0, so that 
/3 > 0. 

The two-dimensional surface of the worldsheet of such a string is spanned by vectors 
ej 1 = X T = (cosh /3, sinh /3, 0, 0) (2.6) 

and 

ef 3) =X CT = (0,0,0,1), (2.7) 
while the vectors 

ej 1 = = (sinh /3, cosh /?,0, 0) (2.8) 

and 

ef 2) =«2 = (0,0,1,0) (2.9) 
are orthogonal to the worldsheet. Vectors (m = 0, 1, 2, 3) form an orthogonal tetrad. 

2. 3. String motion in a weak gravitational field 

The study of scattering of an infinitely long cosmic string by a non-rotating black hole 
is most conveniently carried out in isotropic coordinates, (T, X, Y,Z), for which the 
Schwarzschild metric has the form 

ds 2 = - (1 ~ M/2jR) ' dT 2 + (l + 4;V (dX 2 + dY 2 + dZ 2 ) , (2.10) 
(1 + M/2.R) 2 V 2R J 

where R 2 = X 2 + Y 2 + Z 2 . This metric can be written as 

ds 2 = -(1 - 2$)dT 2 + (1 + 2^)(dX 2 + dY 2 + dZ 2 ) (2.11) 



with 



and 



^ = V^ + ^ 2 + ^ 3 + ^ 4 , (2-13) 
where (y9 = Mj R is the Newtonian potential. 
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The motion of the string in the region with small (f can be written in the form 

X^(T,a)=X^(T,a) + SX"(T,a). (2.14) 

The straight string starts its motion in the plane Y = Y = b. We call it "in-plane" 
(see Figure 1). This plane represents the motion of the free string in flat spacetime. 
Deflection of the string worldsheet from the "in-plane", <5X M , can be decomposed in the 
orthogonal frame as 

= X (m) ef m) • (2.15) 

By linearizing string equations ( |2.3| ) and ( j2.4| ) one gets the following equations for the 

first-order perturbations x @ : 

/ d 2 d 2 \ 

□X(m) = ^-g^ + g^)Xm = fm, (2-16) 

where 

/(o) =2sinh/3cosh 2 /3^ x , (2.17) 
/(i) = 2 sinh 2 (3 cosh (3 (p >x , (2.18) 
/ (2) = -2sinh 2 /5^ y , (2.19) 



and 



/ (3) = -2cosh 2 /3^ iZ . (2.20) 

The solutions of linearized equations ( |2.16| )-( |2.20| ) were obtained and analyzed in 
]. Here we reproduce only a solution for \(2) which describes the deflection of the 
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string in ^-direction, that is in the spatial direction perpendicular to the string motion. 
We assume that a straight string starts its motion with x<2) = at X = — oo then 

Xa(r, a) = -M sinh (3 [H+(t, a) + tf_(r, a)} , (2.21) 

where 

"F 2 + r(r ± a) sinh 2 /?] ^ \{t±<t) tanh/3 



H±(t, a) = arctan 



-arctan 



.(2.22) 



F R sinh /3 

Here i? 2 = r 2 sinh 2 /3 + YJ, 2 + a 2 . This expression for X2(t, cr) can be easily obtained 
from equation (2.22) of Ref. |J if we change r — > r — X / sinh/3 in this expression and 
take the limit X — > — oo. After this redefinition of r, the value of new variable r = 
corresponds to the moment of time when the "distance" R from the string to the black 
hole is minimal. 

The analysis of the late time asymptotic of this solution shows that the string 
approaches the plane Y = Y out , and the shift in position of the string is given by U 

A oo = Y - Y out = 2M 7T sinh (3 = 2M ir 7 v, (2.23) 

where v = tanh/3 is the velocity of the string and 7 = (1 — t> 2 ) -1 / 2 . The late-time 
solution is represented by a kink and anti-kink, propagating in opposite directions at 
the speed of light, and leaving behind them the string in a new "phase" . Furthermore, it 
can be shown that each kink has a characteristic width in the external spacetime given 
by! 

Y 

w = nY cothp = 7i—. (2.24) 

v 

2.4- Scattering data 

The above results were obtained under the assumption that the string always propagates 
in the region where the gravitational field remains weak. In what follows, we describe 
results concerning string scattering in the regime where the gravitational field cannot 

be considered weak. Nevertheless, if the string is not captured it eventually emerges 

I When v — ► second order corrections become important (see [[| 0)- We do not consider the case of 
extremely slow motion of strings in this paper and always assume that v > 0.01c. For a discussion of 
slowly moving strings see [Q. 

Note that there is an additional factor of ir in the expression for the width of the kinks which was 
absent in [pi . 
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from the strong field region and its subsequent evolution is described again by equations 
(|2.16|). The string's late time features are characterized by the amplitude and width 



w of the kinks; we refer to these two parameters as the scattering data. As will be seen 
in the following sections, the breakdown of the weak field approximation manifests itself 
in two ways: the failure of the analytical expressions to correctly reproduce transient 
shapes, which develop when the string reaches periastron, and discrepancies in the late- 
time features of the string, namely the scattering data. 

3. Strong Field Scattering 

3.1. Numerical methods 

In order to study the scattering problem under strong-field conditions (where the weak- 
field approximation breaks down), a numerical solution to the equations of motion must 
be found. To deal with an infinite cosmic string numerically, a finite computational 
domain is required, so a method of truncating the string to a reasonably short segment 
is needed. Such a physical truncation imposes special boundary conditions since it is 
crucial to reproduce the motion of an infinitely long string with sufficient accuracy. 

Two methods of truncating a string were developed. The first places a massive 
particle at each end of the segment of string and requires that the motion of these 
particles mimic the motion of the portions of the string lying outside the region of 
interest. This is accomplished by letting the mass of the end particles go to infinity. It 
should be emphasized that truncating the string and adding infinitely heavy particles 
at its ends is just a technical trick. It allows one to consider finite strings. In making 
the size of the string longer and longer, we finally reach a size where a further increase 
in the string's length does not affect the solution within the chosen accuracy. This may 
be interpreted as saying that the infinitely massive end particles correctly simulate an 
infinite string provided the length of the truncated string is chosen sufficiently long (this 
issue is discussed in Reference 0). 

The second method places the truncation points a reasonable distance away from 
the region of interest, subject to the condition that the end points remain in a weak 
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gravitational field where a perturbative solution to the equations of motion is applicable. 
The analytic weak-field solution is used to prescribe the motion of the boundary points 
of the string. Each method has distinct advantages, and both were used extensively. 
Further, the solvers were tested against one another and agreement under identical 
initial conditions served as an additional validation of the numerical schemes. 

With the choice of internal metric Kab — e 2n VAB, the equations of motion have the 

form 

d 2 X» d 2 X» „„ (dXPdX* dXPdX ,r i\ , . 

on \ ~aZ EZ. ET-^T- =°- l 3 ' 1 ) 



Or 2 da 2 m \ dr Or da da 
The discretization of these equations (described in detail in Reference ||) yields a block- 
tridiagonal system with which the solution over the entire spatial (a) grid is obtained by 
algebraic methods. The boundary conditions, for which X^(t, Cj) = Xf, are imposed 
by replacing the first and last entries in the block-tridiagonal system by the required 
expressions. 

For infinitely massive end particles at the boundaries the following equations are 

substituted into the system, 

d 2 X? „ u dXf dX? . . 

the derivation of this boundary condition is given in Reference [f| . For the perturbative 
boundary conditions, the weak-field solutions of Eqn. ( p7T6| ) are substituted into the 
system; the solutions are given in Reference 0. 

The technical details of the numerical work are not important, but a few general 
comments are in order on the practical aspects imposed by the boundary conditions and 
initial data. 

Under both types of boundary conditions, the accuracy of the numerical results 
may be influenced by the length L of the string segment, meaning that if a string 
segment is too short, the motion of the boundary points induces incorrect behaviour. 
The larger the value of L, the smaller the effect, but the larger the value of L, the 
greater the number of grid points required to adequately resolve the string. The choice 
of string length is therefore a compromise between the need to accurately reproduce 
the behaviour of an infinite string and the need to use computing resources effectively. 
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We found that the massive-particle boundary conditions required a string length of 
at least L = 1000r g (r g = 2M, the Schwarzschild radius of the black hole) to give 
reasonable results, as discussed in Reference |§, whereas the perturbative boundary 
conditions required a much shorter string length, L ~ 100 r g , since the weak-field 
solutions represent the motion of the string boundary points with sufficient accuracy 
at these distances. Therefore, the perturbative boundary conditions greatly reduce the 
number of grid points and time steps needed to resolve the motion of the string near 
the black hole, and hence speed up execution dramatically. 

Initial data takes the form of a straight string, Minkowski solution with the 
stipulation that the initial position of the string is sufficiently far away from the black 
hole to minimize the discrepancy between this approximate initial data and the correct 
initial configuration in Schwarzschild spacetime; the constraint equations are used as a 
guide. The constraint equations are also used to check the accuracy of the calculations 
during the time evolution of the string. 

The evaluation of the discretized form of the constraint equations is carried out 
periodically during the numerical solution. Statistics are computed for the constraints 
(average value over the length of the string and standard deviation) at the current time 
step and reported to an output file. Monitoring that the constraints are consistent 
with zero to several significant digits is done by inspecting this file. If the constraints 
are not satisfactorily maintained, that is if the average value grows unacceptably large 
or undergoes sudden changes, the numerical solution is restarted with new parameters 
(e.g. finer grid, finer step size). The results for the constraint calculation are influenced 
significantly by the initial position of the string. A solution for a string placed initially 
closer to the black hole requires fewer time steps, but generates a poorer result in the 
constraint calculation. Typically, the constraints were expected to be consistent with 
zero from four to six significant figures. Tighter tolerances can also be achieved by 
increasing the number of grid points (and significantly increasing solution time), but 
for the purposes of this study, the above values were deemed a reasonable compromise 
between accuracy and speed. 
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The perturbative solution to the equations of motion of an infinitely long 
cosmic string is most conveniently carried out in isotropic coordinates, (T, X,Y, Z). 
Numerically, however, the scattering problem is best studied in Eddington-Finkelstein 
In-going coordinates since the Christoffel symbols in this coordinate system are 
regular everywhere away from the origin. The Christoffel symbols associated with 
this coordinate system, along with their derivatives with respect to the spacetime 
coordinates, were derived and inserted as analytic expressions in the portions of 
the code dealing with the initialization of the tridiagonal matrices. Also, the 
coded forms of the analytic expressions for the initial data and for the perturbative 
boundary conditions were converted to Eddington-Finkelstein coordinates using 
standard coordinate transformations. 

3.2. Scattering at intermediate velocities 

In this section, we compare the numerical solutions of string scattering to those of the 
weak-field approximation at intermediate velocities (0.1 < v/c < 1) for a range of impact 
parameters. Our basis for comparison is the scattering data. 

Figure 2 shows sections of the numerical worldsheet of a string with initial velocity 
0.76c and impact parameter 40 r g . As shown in Table 1, the maximum amplitude and 
width of the pulses is in good agreement with the weak- field Eqns. (|2.23|) and Q2.24j ). 

Table 1. Scattering data - Schwarzschild weak-field. 



numerical 


perturbation 


Aoo 3.4 r g 


3.6 r g 


w 140r g 


165r ff 


(v = 0.76 c, 


b = 40 Tg) 



When the impact parameter approaches the critical value for capture, the string 
is in a regime where only numerical solutions are available. Figure 3 shows sections of 
the worldsheet of a slow string (v = 0.2c) with near-critical impact parameter (b = 2.5 
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r g compared to Capture = 2.1 r g ). As shown in Table 2, the maximum amplitude of the 
pulses is considerably greater than predicted by weak-field Eqn. ( |2.23p , while there is 
good agreement in terms of the measured width of the kink and Eqn. Q2.24Q . 

Table 2. Scattering data - Schwarzschild strong-field. 



numerical 


perturbation 


Aoo 1.8 r g 


0.95 r g 


w 28 r g 


27 r g 



(v = 0.29 c, b = 2.5r g ) 




Figure 3. Time sequence of string scattering in strong- field regime (numerical results). 
Black hole lies at origin of coordinate system. Initial velocity, 0.29c, impact parameter 
2.5 r B . 



It is important to note that, although the behaviour of the string is qualitatively 
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1.75 




5 10 15 20 25 30 

b/r g 

Figure 4. Comparison of numerical scattering date with weak-field approximation. 



consistent with the weak-field solutions, the total deflection of the string is greater than 
predicted by Eqn. ( |2.23|) . The discrepancy is made clear by plotting the total deflection 



of the string obtained from the numerical solver and comparing it to the prediction 



of Eqn. (|2.23j) . This is done in Figure 4, where the ratio A numerica i / ' A weak is plotted 
for four velocities. It can be seen that numerical and perturbative results converge for 
large impact parameters, and that the transition from weak-field to strong-field occurs 
for impact parameters on the order of 10 r g . Conversely, the curves make it clear that 
the weak-field solutions are quite acceptable down to very small impact parameters. 
The curve which suffers a downturn is the lower velocity one, meaning that Eqn. (|2.23j) 
predicts a greater shift than the numerical data. The velocity of the string in this case 
is at the upper threshold where second-order corrections to the perturbative expansion 
become important. Recall that the expression (|2.23|) was derived using results from 



first-order perturbative calculations. 

4. Ultra-relativistic Limit and Coil Formation 

4-1. Analytical results for ultra-relativistic scattering 

In this section we consider scattering of a straight cosmic string by a black hole in the 
limit where the initial velocity is close to the velocity of light. Numerical studies of the 
evolution of stochastic networks of cosmic strings from initial formation to the matter 
era [|1(| suggest that astrophysical cosmic strings are expected to move with an average 
velocity v ~ 0.7c . For this reason, the results discussed in this section concerning 
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motion at the ultra-relativistic limit are mainly of theoretical interest. Nevertheless, they 
provide a number of analytic expressions against which to validate numerical solutions 
in this regime. 

In the limit v — > c, the gravitational field of a black hole in the reference frame of 
the string takes the form of a shock wave. The corresponding Aichelburg-Sexl metric 
13| is of the form (see Appendix) 



ds 2 = -dX_ dX + + dY 2 + dZ 2 - A^MF5{X + )dX 



+ ■ 



(4.1) 



where F = hip 2 , and 7 = (1 — v 2 ) x l 2 . In these null coordinates, the initial motion of 
the straight string is given by 

(4.2) 



X +\r<0 = T : 



It<0 



t - 2 7 X 



Y \t<0 - Y > 



(4.3) 
(4.4) 



and 



Z\ , n = cr . 

lr<0 



(4.5) 



The problem of scattering of the string in the metric (|4.1| ) was studied earlier (see e.g. 
1^1). For our case, the corresponding solutions are (details in the Appendix) 



and 



X +Uo = T » 



It>0 



T - 27Xo - 27M 



In 1 + 



r + cr 

Yn 



In 1 



8 7 2 M 2 



r + a \ I t — cr 

arctan I —rz — I + arctan 



Y 



r>0 



It>0 



Yq - 2 7 M 



a - 7M In 



Yn 



t + a 

arctan | — I + arctan 



Yn 



Y 

t — a 



Y 







Yq 2 + (r + a) 2 
Y 2 + (t- a) 2 



t — a 



(4.6) 
2 



Yn 



(4.7) 
(4.8) 

(4.9) 
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These results allow one to easily determine the amplitude of the kinks for the ultra- 
relativistic solutions. The maximum amplitude comes from the asymptotic value of 
Y\ 

1 lf>0> 

4»:=y - lim Y\ f>0 = 2n 1 M. (4.10) 

r^oo 

This result is in agreement with the result of weak field approximation ( p.23j ), in the 
limit v — > c. 

4 . 2. Coil formation 

The formation of coils characterises ultra-relativistic scattering. A coil-like configuration 
appears when Z(t, a) is no longer a monotonic function of a for a fixed value of r. This 
occurs when Z CT < 0. From Eqn. ( |4.9| ), one gets 

C+ C- 



2 7 M 



(4.11) 



where (± = r ± a. The function Z t<T has extrema at \(±\ = \Yq\. The minimum of 
Z )(T occurs at ( + = (- = Y , that is at r = Y and a = 0. This minimal value 
Z a = 1 — 2^/M/Yq becomes negative when 

2jM>Y . (4.12) 

This is the condition for coil formation. This condition can be overlaid on the capture 
curve m, as shown in Figure 5. The coil formation region lies above the capture 
curve and to the right of the curve 7 = Y /r g ; it can be seen that coil formation is a 
relativistic phenomenon, and the ultra-relativistic solution predicts that this effect cuts 
off for velocities below v ~ 0.9c 

The maximum width of the coils, w co u, occurs when the solution for Z\~ >0 reaches 
an extremum in both r and a. The maximum width represents the greatest distance, 
relative to the Z = plane, reached by the segment of string forming the coil before the 
coil begins to "collapse" and unwind (see Figure 6). It is straightforward to compute 
the derivatives from Eqn. ( |4.9| ) and show that 



»',„,/ - '2- ■{ r] - ^ln 



1 +y 



W ; ^/!-(^) • (4-13) 



Gravitational Scattering of Cosmic Strings 
10.0 



16 



6.3 



> 4 -° 

Q 

2.5 



Numerical boundary 
of coil region 



Capture curve 



1.6 

y=b/2M 

0.83 0.87 0.91 0.95 

v/c 

Figure 5. Coil formation region. 




Coil region 



4-3. String self-intersection 

The formation of coils is accompanied by string self-intersection. Self-intersection occurs 
if there exists a pair of points on the string worldsheet, (tl, <Ti) and (r 2 , cr 2 ), such that 

X"(Ti,a 1 )=X' t (T2,<T2). (4.14) 

Since this condition forms a system of 4 equations for 4 variables, in the general case the 
self-intersection points are isolated. For the highly symmetric initial-value problem used 
here, a point of self-intersection is located at Z = 0. Really, equation (|4.14j) for X + given 
by ( |4.6|) implies that T\ = t%. Since Y(t, a) given by ( (4.8|) is symmetric with respect to 
cr and as a function of a has minimum at a = 0, the equation Y(t, o~i) = Y(t, a 2 ) has 
solutions o~2 = ±01. Thus for a non-trivial solution of ( 4. 14] ) one must have 

(t 2 ,0 2 ) = (ti,-0i). (4.15) 

For defmiteness, we assume that o~\ > 0. Since Z(t, a) given by ( (4.9|) is an antisymmetric 
function of cr, the function Z vanishes at points of the string self-intersection. 

To find a condition when the self-intersection of the string might occur let us 
consider a function 



F(C) = C - 2 7 Mln [1 + (C 2 /Yo 2 )] . (4.16) 

The point Z = is a point of self-intersection if the equation F{t\ — o~i) = F{r% + o~i) 
has a non-trivial solution with o~\ > 0. This occurs only when the function F is non- 
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monotonic, and the equation 

F'(C) = 1 - 4 7 MC/(F 2 + C 2 ) = (4.17) 

has a solution. The latter condition implies 2jM > Yq. This is exactly the same 
condition of coil formation which was obtained earlier, Eqn. ( |4.12| ). 

Suppose that 27M > Y and denote C± = 27M ± sJa^M) 2 - F 2 . Then the 
function F(() has a local minimum at ( = ( + and local maximum at ( = For any 
value of Fq between -F(C+) and -F(C-) there exist 3 solutions of the equation F(() = Fq 
which we denote by Ci < C2 < C3- Correspondingly, there exist three different branches 
of the solutions of the self-intersection conditions 

n-oi = Ci, n + o-! = ( 2 , 

r 1 -a 1 =( 1 , Tl + a 1 = ( 3 , 

n — ai — C 2 , n + o- 1 = ( 3 . 

It is easy to verify that solutions from different branches have different value of 
T\. Transition between different branches take place when either Ci = C2 = C- or 
C2 = C3 — C+- Simple analysis shows that the self- intersection is stable with respect to 
small perturbations of the initial configuration of the string (before its scattering by the 
black hole). 

Coils are created and disappear at cusps. To shows this let us study the properties 
of vectors X^ T and which are tangent to the worldsheet of the string. Note that 
(X^) 2 = Y\ + Z\ > and hence vector X^ a is spacelike except for points where 
it vanishes. It is easy to see that this occurs either at the point of coil formation, 
(a = 0, r = t_), or at the point of its disappearance, (a = 0, r = r + ), where 

r± = 2 7 M ± ^4(7 M) 2 — Yq . (4.18) 

Constraint equations ( |2.4j) imply 

g» v X^X» = , (X^) 2 = -(X») 2 . (4.19) 

The second relation shows that X^ T is everywhere timelike except for two points where 
X^ a vanishes. It is easy to verify that at these points X^ T is null. Hence, formation of a 
coil and its disappearance always occur at cusps. 
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Figure 6. Time sequence of string scattering in ultra-relativistic regime (numerical 
results). Black hole lies at origin of coordinate system. Initial velocity 0.995c, impact 
parameter 4.0 r g . 

4-4- Numerical results 

Figure 6 shows sections of the numerical worldsheet of a string propagating with an 
initial velocity of 0.995c (7 = 10) and impact parameter of Ar g . The view looking 
down onto the XZ plane show a complicated early phase associated with the evolution 
of short-lived coils and, at late times, two kink-like pulses propagating outward. The 
speed of the pulses is again that of light. The projections onto the YZ plane shows more 
clearly the early evolution of the coils, which indicate that the points on the string near 
the Z = plane undergo a short-lived deflection across the Z = plane. The formation 
of coils suggests that the behaviour of the string at periastron is particle-like in that 
there is insufficient time for tension to play a role, and the points on the string undergo 
a temporary Keplerian deflection. However, once clear of the black hole, the tension can 
again assert itself and solitonic pulses with an S-shaped profile emerge and propagate 
outward. Figure 6 also shows that the string self-intersects in the Z = plane. As the 
coil grows, the intersection point moves gradually downward in the Y-direction. Once 
the coil has reached its maximum size, the S-shaped kinks have formed and begin to 
propagate outward, away from the Z = plane. The intersection point disappears when 
the kinks have completely separated. 

Figure 7 shows YZ sections generated from the analytic weak-field and ultra- 
relativistic solutions for the same physical parameters and intervals of proper time 
as in Figure 6. The two analytic solutions are virtually indistinguishable, the only 
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difference being in the first slice, taken at the time where the string reaches periastron. 
The ultra-relativistic solution shows a completely straight string, as expected, since the 
black hole does not begin to distort the string until it has passed over it. The weak- field 
solution shows a slightly bent string, indicating that even at 0.995c, the influence of 
the black hole is felt before closest approach. For later times, the two solutions are 
indist inguishable . 



perturbation solution 
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Figure 7. Time sequence of string scattering in ultra-relativistic regime (analytic). 
Black hole lies at origin of coordinate system. Initial velocity, 0.995c, impact parameter 
4.0 r a . 



As shown in Table 3, the maximum amplitude of the deflection and the size of 



the coils in the numerical data is slightly larger than predicted by Eqn. ( 4.14 ). (The 
width of coils is difficult to estimate from numerical data, as is consequently shown with 
less precision that the other results.) It is important to note that 7 = 10 for these 
figures, which is low for the ultra-relativistic approximation, but represents a practical 
upper bound for numerical computations^]. Furthermore, the impact parameter is very 
small, and, as shown in the previous section, the approximate solutions are no longer 
completely accurate under these circumstances. 

Since the weak-field and ultra-relativistic approximations are accurate only for 

+ This is due to numerical anomalies that develop when larger 7-factors are considered. These numerical 
anomalies manifest themselves as a rigid deflection of the entire string as it crosses the X — plane. 
The cause of these anomalies was traced to a loss of resolution in the first-order angular derivatives 
d^X^ at large 7-factors. The problem is suppressed by increasing the number of grid points; the large 
grid sizes required result in a numerically intensive solver. 
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Table 3. Scattering data - Schwarzschild ultra-relativistic. 





numerical 


perturbation 




39.7r g 


31.6 r g 


Wcoil 




13.2 r g 




(v = 0.995 c, 


b = 4.00 r g ) 



impact parameters greater than ~ 10r 9 , the boundary of the coil formation is expected 
to shift due to strong field effects. The solid line shown in Figure 5, based on numerical 
tests to detect coil formation, indicates that the boundary is shifted towards lower 
velocities; the strong field near the black hole tends to enhance coil formation^. 

The black hole drives the self-intersection process as its gravitational influence 
draws points on either side of the string towards the Z = plane. What happens 
with the cosmic string forming a cusp which starts growing into a closed coil cannot be 
answered in a simplified model we use. One of the possibility is that after formation a 
small coil, as the result of the string interconnection [|1(J,|T5| the coil creates a loop. In 



this scenario, the black hole driven self-intersection could leave a line of small loops in 
the wake of the string. Another option is that loops never form. Instead blobs of Higgs 
particles are ejected continuously as the budding loop cuts from the string without fully 
forming. This would result in a trail of fastly decaying Higgs particles behind this cusp.| 

5. Near-critical Scattering 

As we already mentioned, there are two possible outcomes when a cosmic string interacts 
with a black hole: either the string is scattered, or it is captured by the black hole. For a 
given initial velocity v the outcome depends on the value of the impact parameter b. In 

other words, in the space of initial data (v , b) there exists a critical line b = b crit (v) which 

* Our numerical results are inconsistent with the findings of Lonsdale and Moss|^|, who observe coils 
at much lower velocities and larger impact parameters. The results discussed here clearly show that coil 
formation is a high- velocity phenomenon. Other inconsistencies with [^) will be discussed in subsequent 
sections. 

jt We are grateful to the referee who attracted our attention to this possibility. 



Gravitational Scattering of Cosmic Strings 21 

separates these two different regimes. We discuss the critical impact parameter for string 
capture in Section 7. This section discusses string scattering where the impact parameter 
is extremely close to the critical impact parameter for capture and the string comes 
close to the black hole. In a similar situation for test particle scattering, the particle 
can execute multiple orbits around the black hole before escaping. Some similarities 
with particle scattering are observed, but there are also important differences that shed 
light on the process of string capture. 

Figure 8 compares a slice of the string worldsheet through the Z = plane 
to the motion of a test particle. A comparison is made of a series of string and 
particle trajectories with identical velocities (v = 0.987c) and nearly identical impact 
parameters. For the particle trajectories, these range from b = 2.60 r g through 
b = 2.65 r g . For the string, the impact parameters are slightly smaller, ranging from 
b = 2.55 r g through b = 2.56 r g . The impact parameters at the lower end of each range 
result in capture, whereas those at the upper end result in scattering. It is easy to see 
that, even at ultra-relativistic velocities, tension influences string motion. Although the 
critical impact parameter for the string is very close to that for a particle, suggesting that 
tension plays a limited role, the dynamics of the string are, nevertheless, still governed by 
tension. This is manifested in the partial winding of the string while gravity dominates 
the dynamics, followed by unwinding of the string as a sufficient interval has elapsed 
for tension to again dominate; the transition between the two regimes is marked by the 
cusp-like feature in box B. By comparison, the escaping trajectories of the test particle 
are far less complicated, showing no folds or cusps. 

The detailed motion of the string in the critical regime close to capture is highly 
sensitive to the initial impact parameter. Figure 8 shows that there are two types of 
capture trajectories, trajectories such as curve 1, where the string crosses the horizon 
directly, and curve 2, where the string folds back on itself before crossing the horizon. 
Curve 2 exhibits a coil-like feature, but does not represent a self-intersection since the 
intersecting points have distinct proper times. Escaping trajectories, such as curve 3, 
have a folded structure. There are two critical cases that mark the transition between 
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each of these three generic curves. The transition between curves of type 1 and type 2 is 
marked by a structure that develops an increasingly cusp-like shape (point B in Figure 
8) as the impact parameter approaches the critical value that marks the transition from 
type 1 to type 2 curves. In practical terms, this impact parameter is difficult to obtain 
since it requires a large number of significant digits (curve 1 and the curve with point 
B have impact parameters that differ by one part in 10 5 ). The transition between type 
2 and type 3 curves is marked by a tangent point, where the string trajectory passes 
twice through the same point (point A) at different times (again, no self-intersection). 
Furthermore, the transition between type 2 and type 3 curves is associated with the 
critical impact parameter for capture. 

A full 3D rendering of the worldsheet nearest the black hole can be generated, but 
it is very difficult to interpret and is therefore not shown here. It is sufficient to note 
that complicated folds and twists are generated as points on the string cross and recross 
while executing their partial orbit of the black hole. At late times, however, these folds 
and twists have dissipated and all that remains are kinks, as shown in Figure 9. This 
figure shows a detailed view near the Z = plane; the kinks evolve and move outward 
along the string, exiting at the top and bottom of Figure 9 and leaving behind the 
straight portion visible at Y « — 15r g . This Figure reinforces the idea that scattering at 
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late times is completely understood in terms of the perturbative solutions. However, it 
may take a significant amount of time for the simple kink/anti-kink picture to emerge. 
In the case studied here, the distortions due to the close encounter persist until the 
string is about 1000 r g past the black hole. 



10 




Figure 9. String worldsheet at late proper time, v = 0.987c. Coils dissipate and kinks 
emerge. 

The numerical results for moderate 7-factors show that the string worldsheet folds 
back on itself, with no evidence of multiple windings or glory scattering as would be the 
case for near-critical scattering of particles. In fact, it seems as if the transition to a 
"direct capture" trajectory (like curve 1 in Figure 8) always occurs before a full turn 
is achieved. The fact that multiple windings are absent is a clear indication that the 
tension in remote parts of the string eventually asserts itself in all cases accessible to 
the numerical solver. The question as to whether the string can complete more than 
one full turn and "wrap" the black hole for larger 7-factors is an open one. 

6. Role of String Tension - Dust Strings 

Since a cosmic string is an extended object under tension, motion of a string near a 
black hole represents the result of the competing influences of tension and gravity. In 
order to shed light on the role of tension in the dynamics of the string, the motion of the 
string is compared to the motion of an array of test particles initially configured with 
the same position and initial velocity as the cosmic string. To do this, consider a family 
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of N test particles arrayed on a line with initial position X^^tq) and initial velocity 

d T X^i(T ), where % is a position index (an integer between 1 and N) that describes the 

initial location of the test particle on the line. These particles each satisfy the geodesic 

equation and constraint, 

d 2 X v i dX p j dX a j _ 



dr 2 p ° dr dr 



d 



dX»idX v i 



= 0. 



dr dr 

It is easy to show that, in flat spacetime, 

Xf (r) = (cosh (P) r, sinh {(5) r + X , Y , a,) (6.2) 

satisfies these equations and represents a straight dust string moving with velocity 
v = tanh f3 where the discrete o~i mimic the a coordinate of the cosmic string. 

It should be emphasized that in the general case dust string dynamics is quite 
different from the dynamics of the Nambu-Goto string. For example, without tension 
a small bump on a dust string simply goes away, while a similar bump on the Nambu- 
Goto string produces two small bumps which propagate away from one another with 
the speed of light. Nevertheless, we found it instructive to compare the motion of a 
straight Nambu-Goto string in the field of a black hole with the motion of a similar dust 
string since it clearly demonstrates the role of string tension in capture and scattering 
processes. 

In the weak-field and ultra-relativistic cases, the role of tension in the string is 
made apparent by comparing the YZ projections of Figures 2 and 6 against those of the 
dust string in Figure 10. Figure 10 shows that a coil always evolves in the dust string 
worldsheet. This coil is due to the Keplerian nature of the trajectories of each particle 
on the dust string. Particles lying on one side of the Z = plane are deflected across 
this plane since their motion is constrained to an orbital plane passing through the black 
hole. Unlike the cosmic string, the Y-axis deflection of the dust string is unbounded 
(this is not to say that the deflection does not reach an asymptotic angular value). Close 
examination of Figures 2, 6, and 10 reveals that, at the early stages of scattering, the 
size of the coils in the dust and cosmic strings are virtually identical, reinforcing the 
idea that tension takes some time to assert itself. Coil formation is a generic feature of 



Gravitational Scattering of Cosmic Strings 25 

200 
150 
100 
50 



-50 
-100 
-150 

-200 

-80 -60 -40 -20 -40 -20 20 40 

Y/rg Y/rg 

Figure 10. Time sequence of dust string scattering in ultra-relativistic regime. Black 
hole lies at origin of coordinate system. Left: Initial velocity 0.995c, impact parameter 
4.0 r g . Right: Initial velocity 0.76c, impact parameter 40 r g . 

the dust string; this is not the case for the cosmic string, where coil formation is subject 
to the low-velocity cut-off effect discussed above. 

In the strong- field regime (with v < c), a comparison to the dust solution is 
uninformative since the impact parameter is well below the critical value for particles 
and the portion of the dust string near the equatorial plane is captured by the black 
hole. However, it does reinforce the idea that internal tension plays an important role 
in the dynamics of the string, helping the string avoid capture for impact parameters 
well below that of the dust string. 

7. String Capture 

In Reference || we presented a plot of the critical impact parameter for capture as a 
function of velocity, the capture curve for cosmic strings. In that paper, we used the 
massive-particle boundary conditions exclusively and showed that, for string lengths 
greater than 2000r g , the capture curve was definitive for velocities greater than 0.2c. 
We are now in a position to revisit this curve using the perturbative boundary conditions 
in the improved numerical solver, and to merge the numerical findings with the ultra-low 
velocity results of to present a complete capture curve. 

The revised capture curve shown in Figure 11 validates two claims made in the 
previous paper. First, that the capture curve for long strings is indeed definitive 
(provided focusing of end points is negligible). This is so because the revised capture 
curve is identical to that of the L = 2000 r g string for v > 0.2c. Second, that there is 
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Figure 11. Cosmic string capture curve for Schwarzschild black hole - data from 
numerical and perturbative results. 



indeed a minimum impact parameter at intermediate velocities, which occurs at v ~ 0.2c 
and has a value of 1.6 r g . However, it is important to note that the new solver is also 
velocity-limited in that, at sufficiently low velocities, second-order perturbative effects 
become important and invalidate the first-order solutions used to set the boundary 
conditions. Nevertheless, analytic results || |7| are available for this velocity range, and 
the numerical results join smoothly with the analytical result for the critical impact 
parameter for very low velocities 



7T 



'capture 




(7.1) 



2\j2v 

and is shown as a dashed line in Figure 11. This capture curve also is in a good 
agreement with a recent result by Page for the critical impact parameter, 



h r JL 

^capture ^ 




tt / t 64 \ 64 ' 



(7.2) 



This formula reproduces the general shape of the curve shown in Figure 11, without the 
fine structure that is revealed by the numerical approach. The error bars f£j] in Figure 11 
are quite small. This suggests that the capture curve is not a smooth line, but has some 
fine structure to it; this is particularly obvious at lower velocities, but this structure 
is also observed on a smaller scale at higher velocities as well. This could possibly 

be explained by the complicated nature of the trajectory of the string for near-critical 

ff As defined in || , the lower error bar denotes the largest impact parameter resulting in capture and 
the upper error bar denotes the smallest impact parameter resulting in escape for a given initial velocity. 
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impact parameters, and its extreme sensitivity to initial data. However, detailed studies 
(carried out subsequent to the drafting of the present paper) of the capture curve at 
ultra-relativistic velocities reveal that the structure is damped out to a certain extent 
by using finer grids, so it is important to note that some (and perhaps all) of the "jitter" 
in Figure 11 may be numerical in origin. 

It should be emphasized that our numerical results demonstrate that in the limit 
v — > c the critical impact parameter reaches the same value, 3\/3GM/c 2 , as the critical 
impact parameter for capture of ultrarelativistic particles. This result is inconsistent 
with the earlier results of Moss and Lonsdale f| . Moreover, the analytical result (|7. 1| ) 
also is inconsistent with the results of Moss and Lonsdale, who quote a low-velocity 
dependence of v~ l obtained from a numerical fit to their capture curve. 

8. Conclusions 

In this paper, we discussed the gravitational scattering of cosmic strings by 
Schwarzschild black holes. This paper brought together earlier results for the critical 
impact parameter for capture, weak-field analytic solutions, and studied the limitations 
of approximate solutions as a function of velocity and impact parameter by generating 
numerical solutions to the equations of motion. Earlier findings have been corroborated 
by these latest results, and numerical, weak-field, and shockwave results are consistent 
with one another in the regions where they overlap. 

The scattering problem was studied in three regimes. In the ultra-relativistic 
regime, we showed that strings form coils, provided that the conditions Yq < 2jM 
and Y > b captU re are satisfied; this phenomenon cuts off for velocities below ~ 0.9c. 
A detailed comparison of analytic and numerical results showed that the perturbative 
results break down for Y < 10r g ; this is a surprisingly small value and suggests that the 
weak-field approximation is quite acceptable for all cases except near-critical scattering. 
Where the weak-field approximation breaks down, in the strong-field regime, kinks are 
produced (without coils for v < 0.9c, with coils for v > 0.9c), but their amplitude is 
larger than predicted by approximate solutions. 

In the case of near-critical scattering, the string worldsheet evolves highly 
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complicated structures. A detailed investigation of string motion in this case showed 
that there are two characteristic trajectories for capture, the first where the string enters 
directly, the second where the string evolves a small coil before crossing the horizon. 
This complicated behaviour, and its extreme sensitivity to initial data, may help explain 
the structure observed in the capture curve; however, a numerical origin to this structure 
has not been ruled out. 

This paper has also shown the critical role performed by tension in determining the 
dynamics of the string. The absence of multiple string windings around the black hole 
and the complicated folds observed in the string worldsheet at near-critical scattering, 
even at moderate 7-factors, are consequences of internal tension. 
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Appendix - Analytical Results in Shockwave Regime 

In this Appendix we present a solution for an ultra-relativistic straight string scattered 
by a black hole. The key observation which leads to a solution of the equations of motion 
in the ultra-relativistic limit is the following: in the reference frame of the string the 



black hole moves with v ~ 1 and its gravitational field is boosted to the shock wave (16 
As a result, before and after crossing the null surface N representing the black hole, 
the string obeys the free equations in flat spacetime. All the information concerning 
the non-linear interaction with the gravitational field of the black hole can be obtained 
in the form of "jump" conditions on the null surface N. Such solutions were studied 
earlier (see e.g. Reference ||14)|). 



To obtain the metric of an ultra-relativistic black hole one starts with the metric 
( p.lOj ) written in the form 

ds 2 = (l + 2^)d^ + 2($ + ^)dT 2 , (A.l) 
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where 

dsl = -alT 2 + dX 2 + dY 2 + dZ 2 , 
and make the boost transformation 

T = i(T-vX), X = >y(X -vT), Y = Y, Z = Z. 
Let X± = f ± X so that 

X + = -y(l-v){T + X), X. = 1 (l + v)(T-X), 
and, since j 2 = (1 — v 2 ) , 
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(A.2) 



(A.3) 



(A.4) 



(A.5) 



X=l{l + v)[ A 



1 - w 
1 + w 



Using Eqns. ( |A.5| ) the metric ( |A.1| ) takes the form 



ds 2 = (1 + 2*) dX + + rfF 2 + dZ 2 ] 

+ £(1 + v) 2 ($ + ¥) ( dX + + 

2 V l + v 



(A.6) 



In the lowest order $ = \Er 
rewritten using Eqn. ( |A.5 ) as 



<p = M/R, where R = VX 2 + Y 2 + Z 2 can be 



X 4 



1 - v 
l + v 



X + ) + -{Y 2 + Z 2 ) 



7 



so that 



7 M(1 - v 2 ) 



1/2 • 



(A.7) 



(A.8) 



(^Y (X + - (g) X_) z + (1 - + ^ 

The metric ( |A.6| ) takes the form 

ds 2 = - dX. dX + + dY 2 + dZ 2 

+ 2y? (rfF 2 + dZ 2 ) + 7 2 y? ((1 + v) 2 dX\ + (1 - v) 2 rfX!) 



(A.9) 



In the limit v — > 1 and 7M fixed one has [14 

hm 70? = hm — , = 
v-+i r ^X 2 + (1 - v> 2 



- 7 M5(X+)lnp 2 



(A.10) 
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where p 2 = Y 2 + Z 2 . From this it follows that the metric in the Aichelburg-Sexl form 
P"5| is obtained 

ds 2 = -dX_ dX + + dY 2 + dZ 2 - 4yMF5(X + )dXl . (A. 11) 

where F = Inp 2 . The boosted metric in this form represents a gravitational shockwave. 
The following components of the Christoffel symbols do not vanish: 

r; + = A 1 MF5\x + ) , r; y = 2Y Y ++ = ^mf ;Y 5(x + ) , 

T- z = 2Y Z ++ = 4 1 MF, Z 5(X + ) . (A. 12) 



ct.(A.13) 



The in-coming straight string motion (|2.5|) in this limit takes the form 



X +lr<0 = T ■ 



X„\ T<Q = T-2yX , Y\ T<0 = Y 



It<0 



First solve the string equation (|2.3| ) in the metric (|A.11|) for the initial conditions 
( |A.13| ) using the conformal gauge in which v '—hh AB = r/ AB , so that □ = —d 2 + d 2 . 
The equation for X + takes a simple form 



DX 4 







(A.14) 



Hence its solution obeying the proper initial conditions is 

X + = t. (A. 15) 

The equation for Y is 

DY = 2 7 MF y 5(r). (A.16) 

The solution obeying the initial conditions ( |A.13|) is solved using Green's Function 
methods HI and has the form 



Y\ 



T>0 



Y - 2 7 M 



T + a \ I r — o 

arctan | — I + arctan 



Y n 



Similarly, solving equation 

DZ = 2jMF z S(t) 



for Z with the initial conditions QA.13D one gets 

Y 2 + (T+o-r 



lr>0 



a - 7M In 



Y 2 + (t - a) 2 



(A.17) 



(A.18) 



(A.19) 
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Expression for X_ can be easily obtained by using constraint equations ( |2.4| ) which 
for the case under consideration take the form (r > 0) 



Y 2 + r 



,2 



z' 2 = X. 



(A.20) 



2(YY' + ZZ') = X'. 



(A.21) 



Here a dot and prime denote derivatives with respect to r and a. By integrating these 
equations one gets 
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